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Abstract 

A high degree of quantum coherence is a crucial requirement for the imple- 
mentation of quantum logic devices. Solid state nanodevices seem particularly 
promising from the point of view of integrability and flexibility in the design. 
However decoherence is a serious limitation, due to the presence of many 
types low energy excitations in the "internal" environment and of "external" 
sources due to the control circuitery. Here we study both kind of dephasing in 
a special implementation, the charge Josephson qubit, however many of our 
results are applicable to a large class of solid state qubits. This is the case 
of 1/f noise for which we introduce and study a model of an environment 
of bistable fluctuatiors. External sources of noise are analized in terms of 
a suitable harmonic oscillator environment and the explicit mapping on the 
spin boson model is presented. We perform a detailed investigation of various 
computation procedures (single shot measurements, repeated measurements) 
and discuss the problem of the information needed to characterize the effect 
of the environment. For a fluctuator environment with 1/f spectrum memory 
effects turn out to be important. Although in general information beyond the 
power spectrum is needed, in many situations this results in the knowledge 
of only one more microscopic parameter of the environment. This allows to 
determine which degrees of freedom of the environment are effective sources 
of decoherence in each different physical situation considered. 
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I. INTRODUCTION 



Solid state coherent system are at the forefront of present day research because of the 
perception that large scale integration may be combined with new physical properties to 
yield new paradigms for nanoelectronics. A concrete example is quantum computation [1-3] 
with solid state devices, with several proposals [4-8] and few recent experiments [9-11]. 
More generally the ability of controlling the dynamics of a complex quantum system would 
open a wide scenario for both fundamental and applied physics. 

Controlled dynamics may be achieved if first it is possible to prepare (write) and measure 
(read) a set of iV observables {Qi} [2]. This set defines the basis of the so called computational 
states | {qi} ) . If moreover it is possible to tune the Hamiltonian of the system, then the 
dynamics of a generic state 

1^,0 = J2 c 9i-?jv(0 \Qi,---,Qn) (1) 

qi-q N 

may be controlled. In general \ip,t) is a superposition of computational states and moreover 
if C{ qi y(t) cannot be factorized in individual functions c qi (t), then \ip,t) is entangled [1,2]. 
Quantum algorithms use superpositions to produce constructive interference towards the 
correct answer and entanglement for the speed up of information processing. Thus certain 
calculations which are practically impossible on a classical computer could be performed 
if coherence could be preserved. Coherence usually denotes situations when a well defined 
relation between the components C{ qi y(t) of the state Eq.(l) exists. In our case this simply 
means that the nanodevice can be described by the pure state Eq.(l), the relation between 
the cs qi \(t) being the solution of the Schrodinger equation. Loss of coherence is due to the 
fact that the Hilbert space of the device is much larger than the computational space [12,13]. 
The system (defined by the set {Qi}) interacts with the environment (defined by all the other 
observables needed to complete the set) [14]. Even a weakly coupled environment may cause 
decoherence [12], i.e. it may destroy the phase relation between C{ qi y(t). The system should 
be described by a density matrix p(t) rather than a pure state as Eq.(l). 

Considering a larger Hilbert space is needed because the nanodevice is a many-body 
object and {Qi} is only a small set of collective variables. Decoherence from these "internal" 
sources represents a serious limitation due to the presence of many low energy excitations in 
the solid state environment. Clever protocols and technological progress may reduce these 
effects, but even in an idealized situation there is an "external" environment of apparata 
(for preparation, measurement, tuning of the Hamiltonian during time evolution) which 
are themselves quantum systems enlarging the overall Hilbert space. To appreciate the 
importance of the external environment, notice that of course we would be happy with a 
system we can manipulate at will, but to obtain easy tunability we have to open a port to 
the external world and this determines decoherence. 
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FIG. 1. (a) Ideal charge-Josephson qubit formed by a superconducting single electron box; (b) 
Charge Josephson qubit in the presence of the circuit and the substrate, with sources of electro- 
magnetic and offset charge fluctuations; (c) Equivalent circuit for the qubit in the electromagnetic 
environment: the impedance Z[u) is modeled by a suitable infinite LC transmission line. 

In this work we consider the charge-Josephson qubit [7]. It is a superconducting island, 
where the charge Q can be stored, connected to a circuit via a Josephson junction and 
a capacitance C 2 (Fig. la). Adding a charge polarizes the surroundings and costs energy, 
provided by the voltage source V x . Then V x may be used to fix Q to the desired value. 
On the other hand Josephson tunneling mixes the charge states. Under suitable conditions 
(charging energy Eq = e 2 /2(Ci + C2) much larger than the Josephson coupling Ej and 
temperatures ksT <C Ej ) only two charge states are important and the system realizes a 
qubit with Hamiltonian 



e Ej 



e{V x ) = AE C {1 - C 2 V x /e) 



(2) 



where the eigenstates { | n) : n = 0, 1} of a z represent a well defined extra number n of 
Cooper pairs in the island and span the computational space. Superconducting qubits [5,6,8] 
are important because they are the only solid state implementations where coherence in a 
single qubit has been observed in the time domain [9,10] and promising experiments exist for 
two-qubit systems [11]. Problems on decoherence can thus be posed in a realistic perspective. 

We will study models of environments which describe decoherence due to fluctuations of 
the external circuit and to 1/ f noise produced by charges which may be trapped close to 
the device (Fig. lb). These have been recognized to be major sources of errors in charge- 
Josephson qubits, because the computational states { | n ) } are coupled to charges moving 
in the environment 1 . The system plus environment Hamiltonian reads 



H = Hq — -Ea z + He 



(3) 



where H.E describes the environment, coupled to the system via the operator E, which acts 
as an extra contribution to the polarization e. 



We will not discuss here other applications of our work, but we stress that charge noise is 
ubiquitous in the solid state and it is important also for solid state qubits based on the electron 
spin, where proposed implementations of two-qubits gates use the Coulomb interaction. Moreover 
the models we consider can be directly applied to flux noise in flux- Josephson qubits [28,29]. 



3 



Since we are interested on measurements on the qubit, the standard road-map is to 
calculate the reduced density matrix p(t) = TrE{W(t)} where W(t) is the system plus envi- 
ronment density matrix. Equations for p(t) can be written in terms of statistical information 
on the environment, i.e. correlation functions like (£ , (t 1 )£ , (i 2 ) ■ ■ )e and system-environment 
correlations, which makes the problem formidable. In practice we hope that in order to make 
predictions we do not really need such a detailed knowledge of the microscopic parameters 
which define 1~Le- This may happen in two remarkable cases, namely if the environment is 
weakly coupled and fast, and if the environment is modeled by harmonic oscillators. Then 
all the information needed on the environment is encoded in the power spectrum of the 
coupling operator E 

dt - (E(t)E(0) + E(0)E(t)) e luJt (4) 

-oo 2 

Roughly speaking if the environment is weakly coupled and fast on the time scales typical 
of decoherence, then the system is unable to probe it in great detail. If the environment is 
made of harmonic oscillators all its equilibrium statistical properties can be derived from the 
power spectrum 2 . Unfortunately this is not enough to describe all the effects of a solid state 
environment, since low-energy excitations may determine memory effects and in general 
specific gates are sensitive to different details of the environment. 

To conclude this introduction we notice that although decoherence comes from the entan- 
glement of the system with its environment, the reduction of the amplitude of the coherent 
signal in specific experiments may often be studied in less fundamental terms. For instance 
one may think to the environment as producing a classical stochastic field which couples to 
the qubit. In this way spontaneous emission is missed but, apart from that, this may be a 
useful point of view for practical estimates of the effect of solid state environments. 



II. SIMPLE ESTIMATES OF DECOHERENCE 
A. Master Equation 

For a weakly coupled environment the trace on the environmental degrees of freedom 
can be approximately calculated in several ways, in second order in the interaction [16]. A 
standard approach leads to the Markovian Master equation in the basis of the eigenstates 
ofW Q [17] 

Pij i^ij Pij ~\~ ^ ' Ti-ijmn Pmn (5) 

mn 

where u>ij is the difference of the energy eigenvalues. The relaxation tensor lZij mn depends 

> 

on combinations of quantities of the kind f£°dt C^ kl (t) and may be calculated if the Green's 
function of the environment iG > (t) = (E(t)E(0)) E , or its Fourier transform 2S(u;)/(l + 



2 The Caldeira Leggett model has been proposed in this spirit to describe the environment in 
Macroscopic Quantum Tunneling and Coherence in Josephson circuits [15] 
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e _|8a; ), is known 3 . In deriving this equation it has been assumed that the environment is 
always at equilibrium and it is fast, i.e. G > (t) should decay on a time scale r c which is 
the smallest scale in the problem. Often the sum in Eq.(5) can be restricted to the secular 
terms (terms such that Uij = u mn ) and the result for the relaxation rate T R , governing the 
exponential decay of the populations pa towards equilibrium, and for the decoherence rate 
which describes the vanishing of the coherences pij, is readily found 

T R =Um 2 6S(n) ; 1^ = T° + ± T R = l - cos 2 S(0) + \ T R (6) 



where Q — ye 2 + Ej is the bare level splitting and tan 6 = — Ej/e is the angle characterizing 
Hq in the Bloch sphere representation. As we will see the major problems in solid state come 
from the so called adiabatic term r°, which depends on low frequencies of the environment. 
The result Eqs.(6) suggests that an optimal operation point is 9 = n/2 and that dephasing 
only depends on the power spectrum of the environment at the operating frequency Q = Ej. 
A typical protocol is depicted in Fig. 2a. 




t time t time 



FIG. 2. Bloch sphere representation of the typical protocol used to generate a phase shift in a 
qubit, by operating with a slightly noisy gate (a) at the optimal operation point, 6 = tt/2, and (b) 
at = 0, when "pure" dephasing occurs. Notice that for charge- Josephson qubit charge states can 
be easily set and measured, so preparation and redout occur along the z axis. 

B. Pure dephasing 

For 6 = the system is sensitive to low frequencies of the environment. This case, usually 
referred to as "pure dephasing", is special in that the Hamiltonian (3) commutes with a z . 
The charge in the island is conserved and no relaxation occurs. However if we prepare 
the qubit in a superposition of charge states the system will dephase and consequently the 
coherences will decay (a typical protocol is depicted in Fig. 2b). 

If the initial density matrix is factorized, W(0) = we{0) <S> p(0), it is possible to write an 
exact expression for the coherences only in terms of the environment [13,18]. 



3 The correlators C< jM (t) read C^ w (t) = Tr E {w E (0)(i\e iHEt H int e- iHEt \j}{l\H int \k)] and C^ w (t) = 
Tr E {w E (0){i\H int \j) (l\e iHEt H int e~ iHEt \k)} , for a factorized initial density matrix W(0) = w E (0)<8 
p(0). The system-environment interaction term for the Hamiltonian Eq.(3) reads Hi n t = —Ea z /2. 
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r(t) = - In Tr E {w E (0) e^^'e-^ 11 } 



(7) 



where H n = H E - (r?/2) 



C. Model for circuit fluctuations 

In order to apply the above results we now specify a model for the environment which de- 
scribes fluctuations of the external circuit, modeled by an external impedance Z{uj) (Fig. lb). 
On the classical level the effect is accounted for by substituting V x — > V x + X(t) in Eq.(2), 
where X(t) are the classical voltage fluctuations X(t) at the impedance. The phenomeno- 
logical quantum model is obtained by applying the following rules. First we substitute 
V x — * V x + X in Eq.(2). Then we introduce a set of harmonic oscillators to describe the 
impedance and to weight quantum fluctuations of X 

which specifies the Hamiltonian (3). Finally the power spectrum of X is identified by the 
voltage fluctuations at Z{uj) 

„ / n , ^ Z(uj) , 3\uj\ w . , 3\uj\ 

S x (uj) = u> 1le K \ ctgh^i = J w ctgh^ 

1 + iuZ{uj)C e ff 2 2 

where C e ff = CiC 2 /(Ci + C 2 ) and to make contact with the standard notation [14] we 
introduced the spectral density of the environment J{uj). This allows to identify S{uj) = 
(4E c C2/e) 2 Sx(u) and to estimate the rates using Eq.(6). At the optimal point the quality 
factor Ej/Tq = (Ci/C e //) 2 2Rq/(ttR) ~ 10 6 for typical values of the parameters (we took 
Z{io) = R and Rq = h/(4e 2 ) is the superconducting quantum of resistence), which would 
allow a single qubit gate to work perfectly. Adiabatic dephasing would lead to an even larger 
Ej/T^ by a factor Ej/K B T. 



D. 1/f noise 



Numerous experiments have shown that the performance of single electron tunneling 
(SET) devices strongly suffer from fluctuations of background charges (BC) located in the 
vicinity of the junctions [19,20]. Their behavior can be visualized as a random extra polariza- 
tion E(t) which produces voltage fluctuations at the device. Voltage noise can be measured 
in SET transistors, and it has been observed to have the 1/f form up to 1 kHz [19]. Using 
experimental data we may identify S{uo) = IQt^AE^/uj where A is about 10~ 6 [19] and 
Eq.(6) at the optimal point leads to a quality factor Ej/T^ ~ 10 3 . This procedure under- 
lies two questionable assumptions: first, one should extrapolate observed 1/f noise up to 
frequencies Ej ~ 10 GHz, second, one should make the ad hoc assumption that there are 



6 



nonequilibrium impurity charges at such frequencies, much larger than the temperature . 
Apart from this Eqs.(6) badly fails in estimating T° which is proportional to the inverse 
of the small (and unmeasurable) low-frequency cut-off 7 m of 1/ f noise. In order to have a 
more reliable estimate it has been proposed that one may describe the environment as a set 
of harmonic oscillators [26]. In this case the exact expression Eq.(7) may be evaluated [13], 
yielding 

T OS c{t) = / — S(u) = . (8) 

JO 71 UJ A 

For times t <C 1/7m, where 7m is the high-frequency cut-off of the 1/ / noise, T osc (t) is 
approximated by T osc (t) 8AEq In (^-) t 2 . The result is still dependent on 7 m in a relevant 
way: each decade of noise, including noise produced by very slow fluctuators, gives exactly 
the same contribution to dephasing. It has been proposed the recipe to use instead an 
ad hoc effective low-frequency cut-off which gives resonable results, but in order to have a 
clear picture of what is going on one should take more seriously the discrete character of 
charge noise. In particular the main difficulty in treating the 1/ / environment by the above 
standard methods is that the large majority of degrees of freedom are much slower than all 
time scales of the evolution of the system, so they give rise to important memory effects. As 
a consequence different gates are sensitive to different details of the environment. 

III. MODEL FOR l/F NOISE 

In this section we introduce a simple model of an environment which yields 1/ / noise, 
which is a suitable set of bistable fluctuators [18]. We first consider the fluctuators to be 
sources of a classical stochastic process, i.e. each fluctuator switches between two con- 
figurations with total rate 7, and determines an extra polarization of the qubit given by 
^Pi(t), where Pi(t) = ±1. The total extra polarization, E(t) = J2i ^Pii^)^ h as power spec- 

OO rj 

trum S{uj) = Y.ii.Vi/2) 2 Jolt (pi(t)pi(0) - 5p ) e iuJt = Ei v 2 /2(l - p 2 ) 7i /(7 2 + to 2 ), where the 

—00 

overline means average on the stochastic processes. The standard assumption [21] of a distri- 
bution of switching rates P(j) oc I/7 for 7 e [y m , 7m] and zero elsewhere leads to 1/ f noise, 
S(u>) = {7r(l — p 2 ) ridV 2 /(4 In 10)} uj~ 1 for frequencies uj G [7™, 7m] (^d is the number of 
fluctuators per noise decade). Already at this level it is clear that the environment presents 
a large number of slow fluctuators, so memory effects are important. 

We introduce a quantum model [18,23] by describing each fluctuator as a localized im- 
purity level connected to a band [24]. The system plus environment Hamiltonian reads 

n = n Q -\a s Y t Vib\b i ( 9 ) 

1 i i 



4 Measurements on the accuracy of single-electron traps showed indirect effects of 1// noise at 
high frequencies [22]. They were interpreted within a similar "one photon" approach wich requires 
the ad hoc assumption that there are nonequilibrium impurity charges at such frequencies. This 
analysis yields the better value A ~ 10 -8 but is certainly non conclusive. 
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Hi = e ci b\bi + Y.l T ki c h b i + h - c -] + Y. £ ki c h c ki ■ 

k k 

Here Tii describes an isolated BC: the operators hi (ftj) destroy (create) an electron in 
the localized level e C i. This electron may tunnel, with amplitude to a band described 
by the operators c^, c ki and the energies Ski- For simplicity we assume that each local- 
ized level is connected to a distinct band. An important scale is the total switching rate 
7i = 27rA/"(e ci )|Tj| 2 (Af is the density of states of the electronic band, \T ki \ 2 ps \Ti\ 2 ), which 
characterizes the classical relaxation regime of each BC. Finally the coupling with the qubit 
is such that each BC produces a bistable extra bias V{. 

It may be argued that the impurity model introduced here describes impurities in metals, 
rather than in insulating substrates. However our aim is to discuss consequences of the 
discrete nature of the BCs and this is the simplest model embedding this feature. Moreover 
this model has been introduced and used in Ref. [23] to successfully explain experiments on 
charge trapping in systems of small tunnel junctions very similar to the charge-qubit. 

Our aim is to investigate the effect of the BC environment on the dynamics of the qubit. 
The picture which emerges from our analysis is that the contribution of the single BC in 
dephasing the qubit depends on the ratio g,i = fj/7j. As a consequence, we distinguish 
between two different kinds of BC: the ones with v/'y <C 1, which we call weakly coupled and 
the ones in the other regime, which we call strongly coupled. Concerning these latter notice 
that we are interested to a physical situation where the Vi are so small that the energy scale 
associated to the total extra bias produced by the set of BCs is much smaller than f2, so 
strongly coupled charges means small 7$. 

Strongly coupled BCs give rise to memory effects and moreover differences of their statis- 
tical properties from those of an oscillator environment (cumulants higher than the second 
are nonvanishing) may be relevant [21]. As a consequence the results Eqs.(6) may be in- 
applicable. One possibility is to calculate higher orders in the Master equation, but even 
if one assumes that the environment is still modelled by harmonic oscillators, the correc- 
tions for 1/f spectrum are of the same order of the leading terms [25]. Here we follow a 
different strategy, namely we enlarge the "system" which allows to treat more accurately 
the dynamics of the BCs, at all orders in the couplings Vj. In the modified roadmap we 
consider only the continuos electron bands in Eq.(9) as the environment and investigate the 
reduced dynamics of the system composed by the qubit and the BCs. We find that weakly 
coupled charges behaves as a source of gaussian noise, whose effect is fully characterized 
by the power spectrum S(uj). On the other hand, from our quantum mechanical treatment 
it emerges that the decoherence due to strongly coupled charges shows pronunced features 
of their discrete character. In the following we will specialize our analysis to the two cases 
= 0,tt/2. 

IV. QUBIT AT THE OPTIMAL POINT 

If the environment is made of a single BC we can implement the new roadmap by solving 
the Master equation Eq.(5) for a system composed by the qubit and the BC. This method is 
unpractical for a large number N of BCs since the number of equations becomes 2 2N — 1. So 
we study the general problem by using the Heisenberg equations of motion. For the average 
values of the qubit observables (a a ), a = x, y, z we obtain (h = 1) 
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N N 

= H v i ( b i b i a y) 5 = ~J2 v i ( b l b i a x) ; (<r z ) = -Ej (a y ) (10) 

i i 



In the rhs averages of new operators which involve also the localized levels and the bands 
are generated and we can write new equations for them. This procedure could be iterated 
and we would obtain an infinite chain of equations, but instead at this stage we factorize 
high order averages, so we are left with a set of 3(N + 1) equations. In practice we ignore 
the cumulants (b\bib^bj) c and (b\bib^bj<j a ) c for i ^ j and insert the relaxation dynamics for 
the BCs in the approximated terms. This method gives accurate results for general values 
of <7j even if Vi/Ej is not very small, as we checked by comparing with numerical evaluation 
of the reduced density matrix of the qubit with one and two BCs. Results are presented in 
Fig. 3 where the time Fourier transform of (a z (t)), proportional to the average charge on the 
island, is shown. We assumed factorized initial condition for the qubit and the BCs. We 
first consider a set of weakly coupled BCs in the range [10~ 2 , 10] Ej which determine 1/ / 
noise in a frequency interval around the operating frequency. The coupling strengths Vi have 
been generated uniformly with approximately zero average and with magnitudes chosen in 
order to yield the amplitude of typical measured spectra [19,20,22] (extrapolated at GHz 
frequencies) . 




FIG. 3. (a) The Fourier transform a z (u) for a set of weakly coupled BCs plus a single strongly 
coupled BC (solid line). The separate effect of the coupled slow BC alone (go = 8.3, dashed 
line) and of the set of weakly coupled BCs (dotted line), is shown for comparison. In the in- 
set the corresponding power spectra: notice that at uj = Ej the power spectrum of the extra 
charge alone (dashed line) is very small. In all cases the noise level at Ej is fixed to the value 
S(Ej)/Ej f» 3.18 x 10~ 4 . (b) The Fourier transform a z (uo) for a set of weakly coupled BCs plus a 
strongly coupled BC (vo/7o = 61.25) prepared in the ground (dash-dotted line) or in the excited 
state (solid line). 

These BCs are weakly coupled, and determine a dephasing rate which reproduces the 
prediction of Eq.(6) (Fig.3 dotted line T^/Ej « 1.65 x 10" 4 ). Now we add a slower (and 
strongly coupled go = i>o/7o — 8-3) BC, in order to extend the 1/f spectrum to lower 
frequencies. The added BC gives negligible contribution to S(Ej) so according to Eq.(6) it 
should not modify T^. Instead, as shown in Fig. 3a, we find that the strongly coupled BC 
alone determines a dephasing rate comparable with that of the weakly coupled BCs. The 
overall F^, is more than twice the prediction of Eq.(6). 

The above result shows that does not depend only on S(Ej). Moreover we checked 
that information beyond the full S(u>) is needed: we considered sets of charges with different 
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N and Vi which realize the same power spectrum S(u) and we find they yield different values 
of r^,. The decoherence is larger if BCs with g >1 are present in the set. In summary these 
results show that Eq.(6) underestimates the effect of strongly coupled BCs and in particular 
that a finite at the optimal point can be obtained even if the 1/ / spectrum does not 
extend up to frequencies ~ Ej. 

If we further slow down the added BC we find that increases toward values ~ 70, 
the switching rate of the BC. This indicates that the effect of strongly coupled BCs on 
decoherence tends to saturate. We discuss later similar results for the case of pure dephasing, 
where this conclusion can be made sharp, but we stress here that this is a consequence of 
the nonlinearity inherent to the discrete nature of the BC environment. In this regime 
we observe also effects related to the initial preparation of the strongly coupled BC (see 
Fig. 3b). As we will discuss later we may describe different experimental procedures using 
different preparations of the environment. Thus we have here a first example of the fact 
that different measurement protocols should be analized separately as far as dephasing due 
to BCs is concerned. 



V. PURE DEPHASING 

In the absence of the tunneling term Eq.(9) is a model for pure dephasing. The analysis 
can be carried out using the exact result Eq.(7). We notice that for our model Eq.(9) 
Ti, v can be decomposed in a sum of commuting terms each referring to a BC, so if we 
assume factorized w E (0) Eq.(7) also factorizes pm(t) = poi(O) ITjli exp{— i(e — Vj/2)t}fj(t) 
in averages referring to a single BC. Using a real-time path-integral technique, the general 
form of fj(t) in Laplace space is obtained [18] 

A -(A)= A + *,,(A)-iy2ftS (n) 

Jl{ X^ + (v 1 /2f + XK 1J (X)+v 1 /2K 2 , 1 (X)- 

where 5p® = 1 — 2(bjbj) t =o specify the initial conditions for the charges. The kernels Kij(X) 
and K2j(X) are expressed by formal series expressions in the tunneling amplitudes Tj. 

An interesting explicit form of the kernels is obtained in the Non Interacting Blip Ap- 
proximation (NIBA) [14], which amounts to approximate the kernels K 1 j(X) and K 2 j(X) at 
lowest order in the tunneling amplitudes Tj. The result is 

Kij(X) — 7j ; K 2J (X)^-^(l/2 + P/27r(X-ie cj ))-i;(l/2 + P/27r(X + ie cj ))} . 

71 

In order to appreciate the physical meaning of the NIBA result, we notice that it is also 
obtained using the Heisenberg equations of motion, which yield a closed set of equations 
provided that the electronic band is assumed in thermal equilibrium. Thus the NIBA result, 
besides describing the dynamics of the qubit plus BC at all orders in the couplings Vj, is 
valid even if the couplings T with the bands are not small. What is possibly missed are 
details of the dynamics of the bands, which is in total agreement with the phylosophy of the 
modified roadmap outlined in Sec. III. 
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A. Nearly incoherent BC dynamics 



We now want to apply these results to a set of BCs which produce 1/f noise. In order 
to produce a clear physical picture, we will consider only the physically relevant limit where 
the BCs have an incoherent dynamics. This limit can be formally obtained from the NIBA 
result by approximating the kernel K 2 j(\) ~ K 2 j(0), and can be expressed in analytic form 

^1 = e~ l£t ft e^ t/2 \A 3 e"'^ 1 -^ + (1 - Aj) e~^ l+a A (12) 
Poi(O) j=1 I J 

where 



otj = y i - g] - 2i gj tgh ((3e cj /2) ; Aj = — (l + a t - i 5p] g 3 ) . (13) 

In the framework of the NIBA this result is valid if e C i , Vi , 7, <C KbT, but again its validity 
is broader. Indeed Eq.(12) it may be obtained in different ways, for instance by analizing 
the system of qubit and BC by a master equation and taking the limit where e C i is the 
largest scale, or even as the exact result of a model where the coupling operator J2iVib\bi 
is substituted by a classical stochastic process E(t) which is the sum of random telegraph 
processes (see App. B). 

The form of Eqs. (12,13) elucidates the different role of weakly and strongly coupled BCs 
in the decoherence process. We focus on the long time limit jt >> 1. Dephasing due to each 
BC comes from the sum of two exponential terms. For extremely weakly coupled charges, 
gj <C 1, whe have aj ~ 1 and only the first term is important, which describes decay of the 
coherences with a rate l/[4cosh 2 (/?£ C j/2)] Vj/^fj. This is precisely the result of Eq.(6) for 
the adiabatic rate For strongly coupled charges, gj 1, each of the two exponentials in 
Eq.(12) expresses roughly the same decay rate oc 7,-, the switching rate of the individual BC, 
and moreover they come with the same weight. The individual <x,s have large imaginary 
parts so the main effect of strong coupling with the qubit is not decay but rather an energy 
shift. The short time limit, which is more important in actual experiments, will be discussed 
in the next subsection. 

The long time limit allows to draw a simple picture of the effect of a BC: the qubit 
is practically insensitive to very fast BCs (gj 3> 1), which are averaged completely. The 
dephasing effect increases as vf /ji as the BC gets slower but eventually saturates. Indeed 
a very slow BC (gj ^> 1) will dephase only when it switches (effect oc 7,) so for most of the 
time it provides a static extra polarization for the qubit. Based on this picture one could 
conclude that slow charges could be neglected, however it is not a priori clear if this is a 
valid choice for the 1/f spectrum, where the number of slow charges is large due to the 
P(7) oc 7 distribution. Eq.(12) allows to answer to this question, but we defer this point to 
section VI and we first discuss in detail other aspects of the result for a single BC. 



B. Single BC 



Results for a single BC are concentrated in Fig. 4a. We consider Eq.(12) for N = 1 and 
look at the quantity 



r(f) 



In 



P01O) 
Poi(O) 
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In Fig. 4 we compare this quantity with the exact result for an oscillator environment, Eq.(8), 
which is reported as the thick dashed line. We show the effect of a single BC for various initial 
conditions and different values of g. Weakly coupled charges (g = 0.1) give a contribution 
close to r osc (t), whereas deviations are observed in the other cases. In particular BCs such 
that g > 1 show slower dephasing compared to an oscillator environment with the same 
S(cu), as a manifestation of saturation effects. Recurrences at times comparable with 1/v 
are visible in T(t). 

Fig. 4a allows to discuss initial conditions of the BC. For each value of g three cases are 
shown. The thick lines correspons to 5p° = 1 (the dotted lines to 5p° = — 1), i.e. the BC is 
initially in the lower (higher) energy classical configuration. These are typical initial situa- 
tions to be considered in a single shot process. In this case T(t) describes dephasing during 
time evolution. Differences in the behavior of T(t) reflect memory effects. In particular for 
7t <C 1 the approximate behavior is T(t) « v 2 t 2 (l — Spl)/8 — 'yvH 3 (1 + 25p 5p — 35po)/24, 
whereas T osc (t) « v 2 t 2 (l — 5p )/8. Thus if 5p° = ±1, T(t) oc t 3 which accounts for the fact 
that a two level system is stiffer than a set of oscillators. On the other hand if we choose 
5p° = 8p, the equilibrium value, T(t) decays more rapidly, T(t) « T osc (t). These latter choice 
of initial conditions corresponds physically to repeated measurements in which we do not 
control the preparation of the BC so T(t) describes both dephasing during time evolution 
and the blurring of the total signal of several realizations of the time evolution with slightly 
different characteristic frequency, a sort of inhomogeneuos broadening. Notice finally that if 
we use the above equilibrium initial conditions T(t) follows T osc (t) for short times, the two 
curves becoming indistinguishable only if moreover g ^ 1. This observation is important 
for the subsequent analysis of inhomogeneuos broadening (Sec. VIII). 

The above analysis of decoherence due to a single BC has clearly evidenced the different 
qualitative influence on the qubit dynamics of weakly and strongly coupled BCs. A single 
weakly coupled, f/7 <C 1, BC behaves as a source of gaussian noise. Thus decoherence only 
depends on the power spectrum of the fluctuator, and does not show any dependence on 
the initial condition of the BC. On the other side decoherence due a single strongly coupled 
BC, v/j 1, displays saturation effects and dependence on the initial condition. 




FIG. 4. (a) Reduced F(t) due to a BC with different preparations (thin lines) for the indicated 
values of g. Inset: longer time behavour for stable state preparation. The thick dashed line 
represents the oscillator approximation; (b) Saturation effect of slow BCs for a 1/f spectrum. 
Relevant parameters (v = 9.2 x 10 7 Hz, = 1000) give typical experimental measured noise levels 
and reproduce the observed decay of the echo signal [20] in charge Josephson qubits. Couplings Vi 
are distributed with dispersion Av/\v\ = 0.2. T(t) is almost unaffected by strongly coupled charges 
(the label is the number of decades included) . 
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VI. 1/F NOISE IN SINGLE SHOT MEASUREMENTS 



We now want to check how the intuitive picture for an environment made of a single 
fluctuator may apply to a set of charges producing 1/f noise. The distribution of BCs 
involves both weakly and strongly coupled fluctuators, so that no typical time scale is present. 
Moreover since a very large number of slow fluctuators is present, it is not a priori clear how 
saturation effects will manifest. 

In Fig. 4b we show the results for a sample with a number of BCs per decade rid = 1000 
and with Vi distributed with small dispersion around the value |i>| = 9.2 x 10 7 Hz. We choose 
initial conditions 5p® = ±1 randomly distributed on the set of N charges to reproduce equi- 
librium conditions in the ensemble, (1/iV) J2j $P j ~ $P- We checked that for every realization 
of the initial conditions we obtained roughly the same overall polarization and the same de- 
phasing. What we are going to calculate is then the result of single shot measurements, 
i.e. the average signal of several experiments on the time evolution of the qubit where the 
total initial polarization of the environment is recalibrated before each experiment. This 
quantity is essentialy dephasing during time evolution and the corresponding protocol is 
such to minimize the effects of the environment. To make reference to a concrete situation, 
the results we present were calculated with parameters close to the experiments [20]. 

In order to illustrate the different role played by the BCs with gj <C 1 and g 3 ; >> 1, we now 
perform a spectral analysis of the effects of the environment. We consider sets of BCs with 
the same 7m = 10 12 Hz 5 and decreasing j m , all producing 1/f noise with the same amplitude 
A in the corresponding frequency range. Solid lines in Fig. 4b represent T(t) calculated with 
Eq.(12). In this example the dephasing is given by BCs with 7j > 10 7 Hz ~ \v 1/10. The 
main contribution comes from three decades at frequencies around \v\. The overall effect of 
the strongly coupled BCs (jj < H/10) is minimal, despite of their large number, showing 
the saturation effect of low-frequency noise. Dashed lines represent T osc (t), the result of the 
approximation of the environment with a set of harmonic oscillators, Eq.(8). In this case 
low-frequency noise does not saturate, each decade of noise equally contributes to dephasing 
and T osc (t) depends on the low-frequency cut-off of 1/f noise. We notice finally that for 
this single shot protocol T(t) is roughly given by T osc (t) provided we use uj ~ |i>| as a low 
frequency cut-off. 

This observation also explains the fact that our results are not very sensitive to the value 
of rid we choose. Indeed in order to reproduce a given amplitude A we must keep constant 
for each decade J^itf ~ n d,v 2 , meaning that the effective low-frequency cut-off ~ |i>| varies 
as n d . For n d — > oo the low- frequency cut-off goes to zero, as in the result Eq. (8) for the 
oscillator environment. 

Finally we point out that even if we performed our calculation by assigning all the mi- 
croscopic parameters of the environment, a reliable estimate of dephasing for single shot 
measurements turns out to depend on a single additional parameter besides the power spec- 
trum S(u>), namely the average coupling \v\ or equivalently the number of charges producing 
a decade of noise, rid- 



5 This value of is surely too large, but in this section we are interested to the way dephasing 
changes by adding low-frequency noise decade after decade 
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VII. COMPARISON WITH THE OSCILLATOR ENVIRONMENT 



In this section we make a more quantitative comparison of our results with the results 
for an equivalent oscillator environment. It is useful to consider our result from the point ov 
view of the classical stochastic approach of App. B which expresses the off diagonal element 
of the reduced density matrix of the qubit as an average over classical stochastic processes 
generated by a set of random telegraph fluctuators. The result for quantum oscillators Eq. (8) 
is the second cumulant expansion of Eq.(12). Alternatively one may verify that 



r osc(t) = ^Y,9i 



d 2 T(t) 



9i=o 



where 5p® = 5p has been posed in T(t). This result is intuitive in that in the limit i>j — > the 
noise is produced by rid oo fluctuators per decade, and its discrete nature is lost. This 
agrees with the fact that in this limit all BCs are weakly coupled. 

To give an idea of the numbers involved we checked this conclusion by calculating T(t) 
from Eq.(12), using different sets of BCs. The results are shown in Fig. 5a. The power 
spectrum S(u) is identical for all the curves shown, which differ in the choice of rid- The 
gaussian behavior is recovered in the long time limit 7 m t ^> 1 for large enough rid (all the 
BCs are weakly coupled). If in addition we take 5p® = 5p, T(t) approaches T 2 (t) also at 
short times. 



VIII. REPEATED MEASUREMENTS AND INHOMOGENEUOS BROADENING 

Measurements of the state of a charge- Josephson qubit involve the measurement of a sin- 
gle extra Cooper pair in the superconducting island. Different strategies of measurements 
have been proposed and experimental prototypes have been produced, but in practice mea- 
surements are not single-shot. In the simple scheme used in Ref. [9,20] the time evolution 
procedure is repeated ~ 10 5 times and what is measured is the total current due to the 
possible presence of the extra Cooper pair in all the repetitions. The signal is then the 
sum over different possible time evolutions of the BCs with initial conditions which are 
also randomly fluctuating. The average effective precession frequency of the qubit is then 
different in different realizations and the overall signal decays faster than a single shot mea- 
surement would yield. In this respect this additional decay of the signal is analogous to 
inhomogeneuos broadening in NMR, where the signal is collected by many noninteracting 
"qubits" each with its specific environment which determines different average precession 
frequencies. One difference may possibly be that strong correlations may exist between the 
different repetitions, due to the strong correlations in time of 1/ f noise. 



14 




10 Hz 




1.5e-09 



FIG. 5. (a) T(t)/T2(t) for a 1/f spectrum between 7„ 



t (sec) 
2 x 10 7 



and 7m 



2 x 10 9 with 

different numbers of BCs per decade: (a) rid = 10 3 , (b) = 4 x 10 3 , (c) = 8 x 10 3 , (d) and 
(dl) rid = 4 x 10 4 , (e) and (el) = 4 x 10 5 . Full lines corresponds to 5p® = ±1, dashed lines to 
equilibrium initial conditions for the BCs. (b) Different averages over 5p® for 1/f spectrum repro- 
ducing the measured noise level of Ref. [20]:|i>| = 9.2 x 10 6 Hz, rid = 10 5 , 7m = 1Hz, = W 9 Hz. 
Dashed lines correspond to the oscillator approximation with a lower cut-off at u = mm{\v\, l/t m }. 

In order to study the overall signal we have to sum Eq.(12) over different sets of {5p° ja } 
representing initial conditions for each repetition starting at t — t a . Equivalently we may 
average Eq.(12) over a suitable distribution of {Sp®} which depends on time. If BCs are 
assumed to be independent this distribution factorizes and since {5p®} enters linearly the 
coefficient A we have only to evaluate Eq.(12) with 5p°j = 5p^(t m ), where 5p®(t m ) is the 
average of the values of 5pj sampled at regular times t a for the overall time t m . As a rough 
estimate we may let 5p®(t m ) = 5p®(0) = ±1 if jjt m < 1 and 5p®(t m ) = 5p for ^jt m > 1, i.e. 
slow charges do not change roughly initial condition (but still they may dephase during time 
evolution) whereas fast charges completely average during t m . In this way the additional 
scale t m enters the problem. 

We consider only the case of long overall measurement time \v\t m ^> 1 which is pertinent 
to the present day experimental conditions. From the results of Sec. VB and Sec. VII we 
may infer that BCs with 7 < l/t m < \v\, being strongly coupled and with 5p® = ±1, are 
ineffective whereas for the other BCs, being averaged, we may take rW(t) w V®.(t) for small 
enough times and, as a result, T(t) jy tm du S(u)(l — cos out) / (kuj 2 ) . This would proof the 
recipe proposed in Ref. [26]. It is intersting to notice that this result and possible extensions 
to higher order effects [25] in the qubit-environment can be obtained with no reference to 
the quantum nature of the environment, and depend on the classical statistical properties 
of an equivalent random process. 

In Fig. 5b, (dotted lines) we show that indeed dephasing calculated as outlined above, 
is roughly given at short times by the oscillator environment approximation with a lower 
cut-off taken at uj ~ min{|t> |, l/t m } (dashed line). We also show results with a different 
averaging procedure 5p Q j(t m ) = l/t m Jq" 1 dtSp(t) which takes into account the strongly cor- 
related dynamics of 1/f noise (solid lines in Fig.5b,). These correlation do not affect the 
results except possibly for t m ~ |i>|. 
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FIG. 6. (a) Bloch sphere descritpion of the echo protocol, (b) Decay of the echo amplitude for 
a 1// spectrum, with noise level as in Fig. (4b) with 7 m = 1Hz and different jm- The dashed lines 
correspond to the oscillator environment approximation. 



IX. CHARGE ECHO 



Echo-type techniques have been suggested [9,26] and experimentally tested [20] as a 
tool to reduce inhomogeneous broadening due to the low-frequency fluctuators of the 1/ / 
spectrum. In the experiment of Ref. [20] the echo protocol consists of a ir/2 preparation 
pulse, a 7r swap pulse and a 7r/2 measurement pulse. Each pulse is separated by the delay 
time t (see Fig. 6). Since the duration of each pulse is negligible with respect to t we can 
estimate the decay of the coherence only looking at the evolution during the delay times. 
Within the semiclassical approach we obtain [30] 



PoiO) 
Poi(0) 



n{4e-^- Rea ^ 7t + A J 2 e-( 1+Rea ^ 7t (14) 



+A 3 e _( - 1_iImaj - )7t + A 4 e~ < - 1+iImaj - )7t | 



where aij is given by Eq.(13) and information on 5p® is contained in Aj. We are interested to 
r( 2 )(t) = -ln|p i(t)/poi(0)|. Again the expansion of Eq.(14) to the second cumulant gives 
the result for an equivalent environment of quantum oscillators 



r£,(0 = - Hdu s(u) 

IT JO 



'i — cosuty 



Results are shown in Fig. 6 for the same set of BCs of Fig. 4 whose parameters are close 6 
to the experimental situation of Ref. [20]. First of all we checked that the dependence of 
r( 2 )(t) on initial conditions 5p® is extremely weak even if we start from out of equilibrium 
configurations of the set of BCs. Moreover r( 2 )(t) ps r^ c (t). This means that the echo pro- 
cedure actually cancels the effect of strongly coupled charges 7 . Only weakly coupled charges 



6 As it is clear from the results of Sec. VI, for a given measured A we had to infer n^. This we 
made by comparing our results of this section with the echo experiment. This allows only a rough 
estimate of because, as in the other cases, the dependence on is weak. 

7 This conclusion is valid as long as the delay time is short, t\v\ <C 1 
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contribute to r( 2 )(t) and the environment behaves as a suitable set of harmonic oscillators. 
We remark that in the regime of parameters we consider, for given noise amplitude the echo 
signal is strongly dependent on the high frequency cut-off, as studied in Fig. 6. 

X. CONCLUSIONS 

In conclusion we have studied dephasing due to charge fluctuations in solid state qubits. 
We analized and compared two models for the environment, namely an environment of 
harmonic oscillators and an environment made of fluctuators with applications to 1/f noise 
which is probably the most serious limitation for these devices. Even if decoherence comes 
in general from the entanglement of the system with its environment, for solid state devices 
measurements of system-environment correlations are extremely hard and the reduction of 
the amplitude of the coherent signal in specific experiments may often be studied in less 
fundamental terms. We discussed the problem of the information needed to characterize 
the effect of the environment. For environments which are weakly coupled and fast or 
environments made of harmonic oscillators the information needed is entirely contained in 
the power spectrum of the operator which couples the environment to the system. For a 
fluctuator environment with 1/f spectrum memory effects and higher order moments are 
important so additional information is needed. We discussed in detail the remarkable case 
of pure dephasing where the reduction of the amplitude of the coherent signal is given by 
correlation functions of the environment alone (see Eq.(7)) and provided exact results for 
the fluctuator environment. In this case the additional information of the environment 
needed depends on the protocol but often reduces to a single parameter. A new energy 
scale emerges, the average coupling \v\ of the qubit with the BCs, which is the additional 
information needed to discuss single shot experiments (alternatively one should know the 
order of magnitude of rid, the number of BCs per decade of noise). For repeated experiments 
the relevant scale is instead min{|t> |, l/t m } where t m is the overall measurement time. Finally 
echo measurements are sensitive to the high-frequency cut-off 7m of the 1/f spectrum. 

Our results are directly applicable to other implementations of solid state qubits. We 
only mention Josephson flux qubits [27,28] wich suffer from similar 1/f noise, originated from 
trapped vortices. In that case the eigenstates of a z are the flux states of the device. Also 
the the parametric effect of 1/f noise on the coupling energy of a Josephson junction [29] 
can be analized within our model, as long as individual fluctuators do not determine large 
variations of Ej. In this case it may be possible that the same sources generate both charge 
noise and fluctuations of Ej. This can be accounted for in our model by choosing the "noise 
axis" as the z axis. 

In this work we have also briefly discussed the possibility of finding optimal operating 
points for the qubit. This idea has been succesfully implemented in the experiment of 
Ref. [10] and consists in operating with external parameters tuned at the points where the 
energy splittings of the system are less sensitive to fluctuations. However a slight deviation 
from the optimal point determines a strong degradation of the performances. An equivalent 
point of view is that the system should be tuned at the point where adiabatic dephasing 
cancels in lowest order, since the effect of a low frequency environment (in particular the 
1/f environment) is minimized. An explicit example is the recent proposal of implementing 
a communication protocol with Josephson junctions [31]. Low-frequency noise can also be 
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minimized by echo techniques, but the flexibility in the implementation of gates is greatly 
reduced. A possibility is to implement quantum computation using Berry phases [8], where 
the design of gates includes an echo procedure, but a detailed analysis of the effect of 1/ f 
noise is still missing. 

Finally we mention that the sensitivity of coherent devices may be used to investigate high 
frequency noise [32]. In particular an accurate matching between measured inhomogeneous 
broadening, echo signal and relaxation may give reliable information on the actual existence 
of BCs at GHz and on the high-frequency cut-off of the 1/ / spectrum. 



APPENDIX A: LINEAR QUANTUM NOISE 

Fluctuations of the electromagnetic circuit can be modeled by coupling the system to 
an environment of harmonic oscillators [14,33] which mimicks the external impedences (see 
Fig.l). In this Appendix we present a model for the electromagnetic environment and we 
derive an effective Hamiltonian H e ff using no phenomeno logical argument [34]. This has two 
motivations. First, since in principle bare circuit parameters are well defined and tunable, 
we want to know precisely how this reflects on H. Second, in general coupling to the 
environment produces decoherence and energy shifts, which may in principle be large. In 
dissipative quantum mechanics shifts are usually treated either by introducing counterterms 
[33] or by writing H e ff in terms of renormalized quantities [33]. The role of induced shifts, 
which is minor in the devices of Refs. [9,10], may be crucial in various situations (e.g. 
geometric quantum computation [8], dynamics of registers and error correction devices). 

We consider the Cooper pair box [5] of Fig. 1. The external impedance is modeled by a 
suitable LC transmission line and the Lagrangian of the system is 



i=l,2 " a 2Z/ a 

where are voltage drops and the Josephson energy is Vj{4>i) — —Ejcos(2e(f>i/h). The 
environment is fully specified by the elements C a and L a . The circuit is introduced by the 
constraint 0i + 02 + J2 a (fia = V x , which allows to eliminate one variable and to write 

' ( i;;(2 ( (/,- 2 <I» //)//) • L b 



2 



where r] = k 2 2 - Ki^i, Cs = 52iQ, C e = CiC 2 /C s , k 1)2 = Ci )2 /C E , and $ = V x - £ a a . 
We next diagonalize L b and obtain the form 



n 



m a . 2 m a u 2 a 2 

rp nr* 



C e V x ^ ] d a x a 



In order to determine the parameters, notice that L5 describes a series C e — Z circuit, 
w hh YspdpXp = Ysa&a = Vz- By comparing the linear response with the known classical 
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dynamics of Vz we determine the spectral density {uj > 0) 



j\uo) = V =Re r zh/c 



1 + wZ(uj)C c 

Notice that since is quadratic this procedure is an exact way to perform the diagonal- 
ization, due to circuit theory, and is valid after quantization due to the Ehrenfest theorem. 
We then perform a (canonical) transformation x = V ~ and obtain the total Lagrangian 
L = L a + Lf, with 

L a = ~7r{x + k> 2 V x - k 2 d a x a f + Vj(2ex/h) 

1 a 

One can verify that the variable canonically conjugated to x is the charge Q in the island. 
The Hamiltonian corresponding to L reads 



Q n d a /2e 

— + Q k 2 2^ — p a - Ej cosl — ; 



' L2m, 



where p a are conjugated to x a and we used the relation ^ a ^/m a = 1/C e . A further 
canonical transformation of the environment (x a = p a / '(m a u a ), p a = —m a uj a x a ) yields 

Q 2 2e 



where c Q = and we introduce J{oS) = 7i"X) a 5(u; ~ ^ajc 2 ^/ (2m a uj a ) = uj 2 J'(u>), the 

spectral density. 

If we isolate the dc bias V x (t) —V x + 5V x (t) and redefine x a + c a C e V x / (m a uj 2 ) — > x a we 
finally obtain 

Q 2 2e 

H eff = T^TT - K 2KrQ + Vj(—X) + Q^2 C « X « 



Notice that the capacitance Ci (and not Cs) enters the Q 2 term, and if we put c a = we 
do not obtain the Cooper pair box hamiltonian. This is correct because the environment 
represents global fluctuations of the circuit, not only of Z. Notice that a static Q shifts the 
equilibrium points of the oscillators and also produces a Q-dependent shift of the zero of 
their energies. This latter can be reabsorbed in the charging energy if we write the oscillator 
hamiltonian using the shifted values x a + c a n 2 Q / {rn a u 2 ), which produces the extra term 
-4Q 2 E a c 2 J(2m a tu 2 a ) = -Q 2 /{2C e ) and 
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H eff = %r - ^V X Q + Vj{^-) + C e SV x V c a x a 
2C S n a 



which reduces to the non dissipative form by letting c a = 0. This is a convenient starting 
point for a weak coupling analysis also because a static shift in the equilibrium points of the 
oscillators has no effect even if it is large. 



APPENDIX B: DEPHASING DUE TO CLASSICAL STOCHASTIC 

FLUCTUATIONS 

In the semiclassical approach the coupling operator Y^i Vibjbi is substituted by a classical 
stochastic process E(t) which is the sum of random telegraph processes. The system is 
initially in a given superposition of eigenstates of 7{q Eq.(2), say | -0, ) = a \a) + (3 \b) . 
We make now a first assumption, namely we neglect relaxation and we write the time 

evolution using the adiabatic theorem | ■?/>, t ) — a \a,t) + (3 eT l h dt ) | b, t ) , where now 
we use the istantaneuos basis of H a d = Hq — \E(t)a z and the instantaneuos level splitting 

Q(t) = sj[s + E(t)} 2 + Ej. The error introduced by the random field E(t) is mostly due 
to the phase factor and we make a second assumption, namely we neglect the difference 
between the istantaneuos basis and the eigenbasis of 1~Lq. In this case the diagonal elements 
of the density matrix \ip,t) {ip,t\ are constant whereas the off diagonal matrix element is 

pUt)=p b a(0)e- i f» dt ' nit ' ) . (Bl) 

At this stage the standard assumption of the semiclassical approach is that averages over 
the environment correspond to averages over the stochastic process E(t). Thus decoherence 
during time evolution is estimated by 



(B2) 

where the overline denotes an average over different stochastic processes with given initial 
conditions. The effect of inhomogeneuos broadening is included by a further averages on 
initial conditions. 

Let us discuss the validity of this result. If Ej — the two assumption leading to Eq.(Bl) 
are rigorously verified and the result is exact. By working out the average we shall pervent 
to the result Eq.(12). In the general case the first assumption is valid if E(t) is slow on 
the time scale set by Q and the second if E(t) Q. The condition of slow environment is 
equivalently expressed by requiring that S(u) is small at frequencies oo ~ fl In this case 
relaxation in leading order (see Eq.(6)) is negligible. Then Eq.(B2) describes the effect of a 
large part of the degrees of freedom composing the 1/ / environment, if not all of them. 

We may apply Eq.(B2) to repeated measurements as described in Sec. VIII. If one as- 
sumes that for small times we can approximate the environment to a suitable set of harmonic 
oscillators then we may approximate Eq.(B2) to the second cumulant. Moreover, under the 



T{t) = - In 



-i J^dt'0(t') 
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assumption of small E(t) it is easy to derive an expression for the decay of the off diagonal 
matrix element 



s in 4 # /■* r 

Tit) = cos 2 6 T osc (t) + —— / dt'dt" \E 2 (t')E 2 (t") — E 2 

81 1 Jo 1 



This expression has been also derived in Ref. [25] where the oscillator environment has been 
studied in detail and we will not pursue this analysis. We only stress that the correction to 
the term containing T osc (t) depends on the so called second spectrum [21] of the environment 
which is known to differ for gaussian and non gaussian stochastic processes. So if the 
correction is important, i.e. if the system turns out to be sensitive to the second spectrum, 
then the approximation of Eq.(B2) to the second cumulant is inaccurate. 

We turn now to the derivation of Eq.(12) for an environment of classical fluctuators. 
Eq.(B2) has to be averaged over the possible realizations of the stochastic process E(t) with 
given initial conditions 8p®. As already discussed we need to perform the calculation only 
for a single Random Telegraph process switching between the two values E a (t) = —v/2 and 
Eb{t) = v/2 with rates 7+ and 7_. (see Fig. (7a)). 
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FIG. 7. (a) Random Telegraph fluctuates switching between valor ±v/2 with rates 7±; (b) 
Each fluctuator can follow four kinds of paths. If intial and final states are equal, the path has an 
even number of transitions, an odd number of flips occurs if initial and final states are different. 

The average E(t) is related to the average population difference between the two state 
a, b of the BC, 5p(t) = (dp - 6pj e -7 * + 5p, by E(t) = -v/28p(t) , where 7 = 7+ + 7_ is the 

total switching rate, 5p° and 5p = (7- — 7+)/7 are respectively the initial and equilibrium 
populaztion difference. For a stationary process 5p° = Sp the Fourier transform of the 
second cumulant [E(t)E{t + r)] c gives the power spectrum of the polarization fluctuations 
v 2 /2{l-p 2 ) 1 /{ 1 2 + u 2 ). 

In order to calculate the generating functional Z(t) = exp{— % J * dt'E(t')}, we have to 
take into account all the possible realizations of the stochastic processes E(t) with fixed 
initial conditions for the fluctuator, Sp°. To this end we have to consider that the fluctuator 
can follow one of the four kinds of paths illustrated in Fig. 7b. 

Then we have to perform (i) a sum over all possible number of switches, (ii) an inte- 
gration over the corresponding transition times. For each process in Fig. 7b we can define a 
probability density ir at @(t; t m — ti\0),(a, f3 can be either a or b) of finding the fluctuator in the 
state P at time t once it has been prepared in the state a at time t — after m transitions 
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at times ti, i — 1, ...m. The probability densities of the elementary processes of having at 
most a single transition in a fixed time interval are easily found from the statistics of the 
populations of the two states a and b, p a (t), a — a,b. The probability densities to stay 
in the states a and b are 7r aa (t 2 \ti) = ex p[ _ 7+(^2 _ ti)], and iTbb(t 2 \ti) = exp[— 7 _(t 2 — ii)]- 
The probability density of a single flip at time t\ in the interval t — to for instance from 
a to b reads 7i a b(t; ti\t ) = exp[— j-(t — £i)]7+ exp[— j+(ti — to)]. With this, the probability 
density of the complex processes 7r Q)( g(i; t rn — ti|0) are given by: 

n aa (t;t 2N -t^O) = ( 1+1 _) N e-^ T ^- T » = ( 7+7 _)^"7 + *-7-^T 2Jv ^ 
7i ab (t; t 2N+1 - * x |0) = { 1+1 -) N 1+ e-i+ T ^- T » = ( 7+7 _) JV 7+ e^- t -^^ T -+ 1 (B3) 

The probability densities 7Tb&(£; £ m — ti|0), nba(t', t m — ii|0) are found from the expressions 
n a a{t', t m — t\\0) and 7Tb a (t; t m — ti|0) with the replacements 7+ — > 7 _ and 5p — > — 5p. Each 
K a ,j3{t'i~tm ~ 1 1 1 0) only depends on the number m = 2iV, (2A^ + 1) of switches and on the 
total times T a during the random fluctuator has assumed value E a . Thus we have: 

z(t) = J2 P ° a z a At) 

*«,/3(*)=E / dt m- dt 1 ir a ^t;t m -t 1 \0)e-< E ' Ta+E ^\ (B4) 

where p° is the probability to find the fluctuator in the state a at time t — 0, and m = 2A?" 
if a — (5 (the term m = being 1), m = 2A^ + 1 if a ^ (5. The total times spent by the 
fluctuators in the state a or b, T a and TJ,, are given respectively by 

2N 

T b = J2 (-l)X = T 2JV T a = t-T 2N a = (3 

n=l 
2N 

T b = J2(-l) n tn + t-t 2N+1 = t + T 2N+1 T a = -T 2N+l a^P (B5) 

n=l 

From eqs. (B3), (B5) Z(t) can be written in the simple form 

Z{t) = Ae-^ {l ~ a)t + (1 - A)ei {1+a)t (B6) 

where A and a are defined in Eq. (13). An interesting approximation of Eq.(B6) is obtained 
by taking its second cumulant expansion. This amounts to consider a Gaussian process 
E(t) having the same power spectrum S(cu). Estimating the average given in Eq.(B6) by its 
second cumulant and taking 5p° = 5p we get Eq.(8). 

The above exact results of the semiclassical analysis represent an important limiting 
case of decoherence due to a quantum discrete environment. Eq. (B6) reproduces in fact 
the results presented in Section VA for the pure dephasing decay due to an environment of 
quantum bistable fluctuators, in the regime in which each fluctuator performs a relaxation 
dynamics. 

We thank L. Faoro, A. D'Arrigo and F. Taddei who collaborated to part of the work 
presented here. We warmly thank G. Giaquinta who initiated the work on mesoscopic physics 



22 



at the University of Catania. We thank M. Palma, A. Shnirman, G. Schon and C. Urbina 
for discussions which greatly sharpened the point of view presented in this work. Very 
useful discussion with L. Amico, J. Clarke, D. Esteve, F. Hekking, P. Lafarge, J. Friedman, 
M. Grifoni, D. van Harlingen, J.E. Mooji, Y. Nakamura, Y. Nazarov, F. Plastina, J. M. 
Raimond, V. Tognetti, U. Weiss, D. Vion and A. Zorin are finally acknowledged. 



23 



REFERENCES 



[1] A. Ekert and A. Jozsa, Rev. Mod. Phys., 68, 733 (1996); Quantum Computa- 
tion and Quantum Information Theory, edited by C. Macchiavello, G.M. Palma, A. 
Zeilinger, World Scientific (2000). 

[2] M.Nielsen and I.Chuang Quantum Computation and Quantum Communication, Cam- 
bridge University Press, (2000). 

[3] Experimental implementation of quantum computation, edited by R.G. Clark, Rinton 
Press, Princeton (2001). 

[4] D. Loss and P. Di Vincenzo Phys. Rev. A, 57, 120 (1998). 

[5] Y. Makhlin, G. Schdn and A. Shnirman, Rev. Mod. Phys., 73, 357 (2001). 

[6] D.A. Averin Sol. State Comm., 105, 659 (1998); L.B. Ioffe et al. Nature, 398, 679 
(1999); J.E. Mooij et al. Science, 285, 1036 (1999). 

[7] Y. Makhlin, G. Schdn and A. Shnirman Nature 398, 305 (1999); A. Shnirman, G. Schdn 
and Z. Hermon Phys. Rev. Lett. 79, 2371 (1997). 

[8] G. Falci et al. Nature, 407, 355 (2000). 

[9] Y. Nakamura, Yu.A. Pashkin, J.S. Tsai Nature, 398, 786 (1999). 
[10] D. Vion et al. Science296, 886 (2002); Y. Yu et al. Science, 296, 889 (2002); J. Martinis 
et al. Phys. Rev. Lett., 89, 117901 (2002); J. Friedman et al. Nature, 406, 43 (2000); I. 
Chiorescu et al. private communication. 
[11] Yu. A. Pashkin et al. cond-mat/02 12314. 
[12] W. Zurek Physics Today, 44, 36 (1991). 

[13] G.M.Palma, K.-A.Suominen and A.K.Ekert Proc. Roy. Soc. London A, 452, 567 (1996). 
[14] U. Weiss Quantum Dissipative Systems 2nd Ed (World Scientific, Singapore 1999). 
[15] A. Leggett et. al Rev. Mod. Phys., 59, 1 (1987). 

[16] A. G. Redfield IBM J. Research Develop, 1, 19 (1957); M. Lax Phys. Rev., 145, 110 
(1966). 

[17] C. Cohen- Tannoudji, J. Dupont-Roc and G. Grynberg Atom-Photon Interactions, 

Wiley-Interscience (1993) 
[18] E. Paladino et. al. Phys. Rev. Lett., 88, 228304 (2002). 
[19] A.B. Zorin et. al Phys. Rev. B, 53, 13682 (1996). 
[20] Y. Nakamura et. al. Phys. Rev. Lett., 88, 047901 (2002). 
[21] M.B. Weissman Rev. Mod. Phys., 60, 537 (1988). 
[22] M. Covington et al. Phys. Rev. Lett., 84, 5192 (2000). 
[23] R. Bauernschmitt and Y.V. Nazarov Phys. Rev. B, 47, 9997 (1993). 
[24] G.D. Mahan Many-Particle Physics Kluwer Academic, New York, (2000) 
[25] A. Shnirman,Y. Makhlin, G. Schdn, Physica Scripta, T102, 147, (2002). 
[26] A. Cottet et al. in Macroscopic Quantum Coherence and Quantum Computing edited 

by D.V. Averin, B. Ruggiero and P. Silvestrini, (Kluwer Pub., 2001),pg.lll. 
[27] L. Tian et al. in Proceedings of the NATO-ASI on Quantum Mesoscopic Phenomena 

and Mesoscopic Devices in Microelectronics, edited by I.O. Kulik and R. Elliatioglu 

(Kluwer Pub. 2000), pg. 429. 
[28] J.E. Mooij et al, Science, 285, 1036 (1999). 

[29] D. J. Van Harlingen et al. in Quantum Computing and Quantum Bits in Mesoscopic Sys- 
tems Proceedings of the International Workshop on "Macroscopic Quantum Coherence 
and Computing", Napoli 3-7 Giugno 2002 (Kluwer Academic Plenum Press) (2003). 



24 



[30] E. Paladino et al. in Quantum Computing and Quantum Bits in Mesoscopic Systems, 
Proceedings of the International Workshop on "Macroscopic Quantum Coherence and 
Computing", Napoli 3-7 Giugno 2002 (Kluwer Academic Plenum Press) (2003). 

[31] F. Plastina and G. Falci, Phys. Rev. B, 67, 224514 (2003) 

[32] R. Aguado and L. P. Kouwenhowen Phys. Rev. Lett. 84, 1986 (2000) 

[33] A. O.Caldeira and A. J. Leggett Ann. Phys.149, 374 (1983). 

[34] E. Paladino et. al Physica E 18, 39 (2003) 



25 



